The spherical 2 + p spin glass model: an analytically solvable model with a 

glass-to-glass transition. 



in 
o 
o 

(N 



o 

CD 



A. CrisantJ3 and L. Lcuzzjj] 

Dipartimento di Fisica, Universita di Roma "La Sapienza" 
Istituto Nazionale Fisica della Materia, Unita di Roma, 
and SMC, P.le Aldo Mow 2, F00185 Roma, Italy and 
Istituto Studi della Complessitd (ISC), CNR, Via dei Taurini 19, 1-00185 Roma, Italy 

We present the detailed analysis of the spherical s + p spin glass model with two competing in- 
teractions: among p spins and among s spins. The most interesting case is the 2 + p model with 
p > 4 for which a very rich phase diagram occurs, including, next to the paramagnetic and the 
glassy phase represented by the one step replica symmetry breaking ansatz typical of the spherical 
p-spin model, other two amorphous phases. Transitions between two contiguous phases can also be 
of different kind. The model can thus serve as mean-field representation of amorphous-amorphous 
transitions (or transitions between undercooled liquids of different structure). The model is analyt- 
ically solvable everywhere in the phase space, even in the limit where the infinite replica symmetry 
breaking ansatz is required to yield a thermodynamically stable phase. 

PACS numbers: 75.10.Nr, ll.30.Pb, 05.50.+q 
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Spin glasses have become in the last thirty years the 
source of ideas and techniques now representing a valu- 
able theoretical background for "complex systems" , with 
applications not only to the physics of amorphous mate- 
rials, but also to optimization and assignment problems 
in computer science, to biology, ethology, economy and 
finance. These systems are characterized by a strong de- 
pendence from the details, such that their behavior can- 
not be rebuilt starting from the analysis of a 'cell' con- 
stituent but an approach involving the collective behav- 
ior of the whole system becomes necessary. One of the 
feature usually expressed is the existence of a large num- 
ber of stable and metastable states, or, in other words, a 
large choice in the possible realizations of the system and 
a rather difficult (and therefore slow) evolution through 
many, details-dependent, intermediate steps, hunting its 
equilibrium state or optimal solution. 

Mean-field models have largely helped in comprehend- 
ing many of the mechanisms yielding such complicated 
structure and also have produced new theories (or com- 
bined among each other old concepts pertaining to other 
fields) such as, e.g., the spontaneous breaking of the 
replica symmetry and the ultrametric structure of states. 

Among mean-field models spherical models are ana- 
lytically solvable even in the most complicated cases. 
Up to now only spherical models with one step Replica 
Symmetry Breaking (1RSB) phases were studied, mainly 
due to their relevance for the fragile glass transitioniS*^ 
The possibility of the existence of Full Replica Symme- 
try Breaking (FRSB) phases in spherical models was first 
pointed out by Nieuwenhuizeni on the basis of the sim- 
ilarity between the replica free energy of some spherical 
models with multi-spin interactions and the relevant part 
of the free energy of the Sherrington-Kirkpatrick (SK) 
modeli^ A complete analysis, however, was not provided 
up to now. The problem has been considered some years 
later— in connection with the possible different scenarios 
for the critical dynamics near the glass transition, 8 there- 



fore analyzing only the dynamical behavior in the 1RSB 
phase. 

The model we present here, the 2 + p spherical spin 
glass model, displays four different phases: together with 
the replica symmetric, the 1RSB and the FRSB phases 
also a new phase occurs. The evidence for the existence of 
such peculiar amorphous phase has been first presented 
in Ref. |9| and, for what concerns the organization of the 
states, it seems to yield the properties of a glass up to 
the first level of the ultrametric tree (i.e. inside and just 
outside a valley of the free energy landscape) and those 
of a spin-glass above. 

Concentrating on the study of amorphous materials, 
in recent years some evidence has been collected for the 
existence of amorphous to amorphous transition (AAT), 
in certain glass-forming substances. One way of looking 
at an AAT has been to consider the kinetics of the coor- 
dination transformation occurring in strong glasses such 
as the vitreous Germania (Ge02 - from fourfold to six- 
fold coordination rising the pressure) and Silica (Si02 - 
from tetrahedral to octahedral coordination)^ Exactly 
as for the liquid-glass transition also this transition is not 
a thermodynamic one, but it amounts to a qualitative 
change of the (slow) relaxation dynamics, apparently ex- 
pressing a recombination of the glass structure (see also 
the numerical simulations of Ref. [ll| for a different point 
of view). Another kind of pressure induced AAT takes 
place in densified porous silicon, where the high-density 
amorphous Si transforms into a low density amorphous 
Si upon decompression^ A similar transition also takes 
place in undercooled water^ 

Theoretical models have been introduced to describe 
an AAT. As, for instance, a model of hard-core repul- 
sive colloidal particle subject to a short-range attrac- 
tive potential that induces the particle to stick to each 
^] ler| i4 i i5 i i6 j n j.jjg f ramewor k of the Mode Coupling 

Theory (MCT) it has been shown that the interplay of 
the attractive and repulsive mechanisms results in the 
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existence of a high(er) temperature "repulsive" glass, 
where the hard-core repulsion is responsible for the freez- 
ing in of many degrees of freedom and the kinetic ar- 
rest, and a low(er) temperature "attractive" glass that 
is energetically more favored than the other one but 
only occurs when the thermal excitation of the parti- 
cles is rather small. Such theoretical and numerical pre- 
dictions seem to have been successfully tested in recent 
experimentsiiiiiLiS Another model where AAT is found 
is the spherical p-spin model on lattice gas of Caiazzo 
et alOU where an off-equilibrium Langevin dynamics is 
considered, thus going beyond the MCT assumption of 
equilibrium. The model we consider here might, as well, 
be a good mean-field representative of an amorphous- 
amorphous transition. 

The goal of this paper is to give a detailed discussion 
of the different solutions describing the low temperature 
phase of the spherical 2 + p spin glass model. As it hap- 
pens in systems with a phase described by a 1RSB solu- 
tion one must distinguish between the "static solution" 
obtained from the partition function and the "dynamic 
solution" obtained from the relaxation dynamics j 2 ^* 2 ^ 
To keep the length of the paper reasonable we shall con- 
sider in detail only the static approach and introduce the 
the dynamic solution with the help of the complexity^ 
The complete dynamic approach will be presented else- 
where. 

In section [I] we present the spherical 2 + p spin-glass 
model. In section [n] its static behavior is studied with 
the help of the replica trick 24 and the Parisi replica sym- 
metry breaking scheme: 25 four different phases occur, to- 
gether with the relative transitions between them. The 
nature of the phases is thoroughly discussed and analyt- 
ical exact solutions for order parameters, transition lines 
and thermodynamic functions are provided all over the 
parameter space. In section IIIII the existence of an ex- 
ponential number of energetically degenerate pure states 
is considered by analyzing the complexity function. The 
connection between the "marginal condition" (maximum 
of the complexity in free energy) and the dynamical so- 
lution leads, in section IIVI to the discussion of the lat- 
ter in those cases where it differs from the static one. 
In Appendix A and B we show, respectively, the Parisi 
anti-parabolic equation for the 2 + p model, and its an- 
alytical solution. In appendix C some basic features of 
the behavior of the much simpler s + p model (s, p > 2) 
are given£2§i Eventually, in appendix D, a proof is given 
that no phases other than those here presented exist for 
the spherical 2 + p spin glass model. 



I. THE MODEL 

The spherical 2 + p spin glass model is defined by the 
Hamiltonian 
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where p is an integer equal or larger than 3 and cr, are N 
continuous real spin variables which range from — oo to 
+oo subject to the global spherical constraint 



A' 



N. 



(2) 



The coupling strengths J^' i (p = 2, 3, . . .) are quenched 
independent identical distributed zero mean Gaussian 
variables of variance 
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The scaling with the system size N ensures an extensive 
free energy and hence a well defined thermodynamic limit 
N — > oo. Without loosing in generality one may take 
either J2 or J p equal to 1 since this only amounts in a 
rescaling of the temperature T. To keep the discussion 
as simple as possible in this paper we shall not consider 
the effect of an external field coupled linearly with the 
spin variables <Ji. 

The properties of the model strongly depend on the 
value of p. For p — 3 the model reduces to the usual 
spherical p-spin spin glass model in a field 2 with a low 
temperature phase described by a 1RSB solution. For 
p > 3 the model exhibits different low temperature 
phases which, depending on the temperature and the ra- 
tio Jp/Ja between the strength of the non-harmonic and 
the harmonic parts of the Hamiltonian, are described by 
1RSB and/or FRSB solutions. 



II. THE STATIC SOLUTION 

The static solution is obtained from the minimum of 
the free-energy functional computed from the partition 
function. The model contains quenched disorder and 
hence the partition function must be computed for fixed 
disorder realizations^LSSiSi 

Z N [J^\ J(p>] = Tr CT exp(-/3H[j( 2 \ J^;a}) (4) 

with (3 = 1/T. We have explicitly shown the dependence 
of the Hamiltonian on the realization of the random cou- 
plings to stress that Zm is itself a function of the cou- 
plings realization. The trace over the spins is defined 
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Tr CT = 2VN 
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(5) 



and includes the spherical constraint As a conse- 

quence Tr CT (l) is equal to the surface of the N dimen- 
sional sphere of radius TV 1 / 2 and its logarithm gives the 
entropy of the model at infinite temperature. 

The partition function Z/v is a random variable, there- 
fore the quenched free energy per spin is given by 



N 



In Z N [JW, J(p)] 



(6) 
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where here, and in the following, (• ■ ■) denotes the average 
over the realizations of all couplings in the Hamiltonian: 



(...)= / dP[jW]dV[J^] (•••) 



(7) 



The thermodynamic limit N — > oo of the free energy 
$ = lirn/v^oo $at is well defined and is equal to the limit 
-limjv^ooln^jv[J (2) , J {p) ]/ N(3 for almost all coupling 
realizations (self-average property) . 

The analytic computation of the quenched free energy 
i.e., of the average of the logarithm of the partition func- 
tion, is quite a difficult problem, even in simple cases as 
nearest neighbor one dimensional models. Since the in- 
teger moments Z^ of the partition function are easier to 
compute, the standard method to evaluate © uses the 
so called "replica trick" by considering the annealed free 
energy <i>(n) of n non-interacting identical 'replicas' of 
the system) 28 i 29 i 30 
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The quenched free energy <E> is then recovered as the con- 
tinuation of <f>(n) down to the unphysical limit n = 0)21 



$ = — lim lim 

N^oo n—tO 



(Z N [jW,J(p)}) n -1 



Nf3n 



lim $(n). (9) 

n— >0 



In the last equality we assumed that the replica limit 
n — > and the thermodynamic limit N — > oo can be ex- 
changed. The existence of such a limit has been recently 
rigorously provedi 32 ' 33 

The replica method gives a simple way of perform- 
ing the disorder average, at the expenses of introducing 
an effective interaction among different replicas in the n- 
dimensional replica space. The interested reader can find 
a detailed presentation of the replica method for disor- 
dered systems in Refs. [2^ and [291 and for the particular 
case of spherical models in Ref. |2j. 

Applying the replica method the integer moments of 
the partition function of the spherical 2 + p spin glass 
model can be written, neglecting all unnecessary con- 
stants and terms irrelevant for N — > oo, as^ 



Zj = e nNs ^ I dq a0 e NG M (10) 



<J>0 



Q</3 



where s(oo) = (1 + ln27r)/2 is the entropy per spin at 
infinite temperature and G[q] the functional: 

1 1,n 1 
G [<?] = ^$>(<M) + ^detq (11) 

(12) 



a0 



g{x) = —x + -^x p . 
2 p 

with fi p = ((3J p ) 2 p/2. We also introduce the two addi- 
tional functions 

A( x ) = ± g ( x ), E(x) = ^A(x) (13) 



whose utility will be clear in a short while. 

The symmetric n x n real matrix q a p is the replica 
overlap matrix 



1 



N 



(14) 



The spherical constraint, Eq. @, implies that the di- 
agonal elements of the matrix q are all equal to one: 
q aa = q = 1- 

In the thermodynamic limit N — > oo the integrals in 
(fTUfl can be evaluated by the saddle point method and 
the quenched free energy per spin $ reads: 



0$ = s(oo) + lim -G[q] 

n— > n 



(15) 



where G[q] must be evaluated on the solution of the sad- 
dle point equation which, in the n — > limit, reads: 



HQap) + (<? l ) a p = 0, a ^ (3. 



(16) 



Stability of the saddle point calculation requires that 
the quadratic form 

- £ E(« a/S ) (<5g^) 2 + Tr^- 1 ,5q) 2 (17) 

a/3 

must be positive definite. Here above Sq a /3 = 5q a p (= 
8q/3a) is the fluctuation of q a p from the saddle point value 
(|16[) . Details of the derivation of these equations can be 
found in Ref. 0. 

The structure of the overlap matrix q a p reflects the or- 
ganization of the different thermodynamic states, called 
pure states, in which each replica can be found. This 
however does not follow from the replica calculation and 
therefore to evaluate explicitly G[q] some ansatz on the 
structure of q must be imposed. 



A. The Replica Symmetric Solution (RS) 

The simplest ansatz is the one in which all replicas are 
in the same pure state, so that g Q( a cannot depend on the 
replica indexes: 



q a p = (1 - q)5ap + q- 



(18) 



This is called the Replica Symmetric (RS) ansatz. This 
assumption is reasonable for coupling strengths not too 
large or high temperatures, i.e., fi2 and /i p small enough. 
In both cases, indeed, the system can explore almost the 
whole available phase space so that different replicas will 
be found in the same pure state. 

Inserting the RS form (|T%|) of q a p into Eq. ijTT)) one 
gets 

2 lim - G[q] - 5 (1) - 9(q) + ln(l - q) + (19) 
n->o n 1 — q 
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where q is the solution of the RS saddle point equation 

A{q) ~ (T^^ 2 ^^ 1 - (T^ = - (20) 

In absence of external fields the saddle point equation 
always admits the "paramagnetic" solution q = 0. How- 
ever since g/(l — q) 2 diverges as q — > 1 and vanishes for 
q = for particularly chosen values of the parameters \x v 
and fi2 there may also be solutions with < q < 1. 

The RS solution is stable, i.e. the quadratic form 
(|17|l is positive definite for n — > 0, provided that the 
eigenvalue^ 



Ai 



■S(g) 



"M2 



1 



(21) 



is positive. If q ^ then dividing the saddle point equa- 
tion l|20(l by q and adding the result to l|21|l one gets that 
the requirement Ai > is equivalent to 



u p (p - 2) qP- 2 > 0. 



(22) 



This inequality cannot be satisfied for any q > 0, thus we 
are left with the q = solution only. 

For q = the eigenvalue Ai reduces to Ai = 1 — //2- 
Therefore, in the (u p ,U2) plane the paramagnetic RS so- 
lution q — is stable everywhere below the \ii — 1 line, 
which represents the De Almeida-Thouless lineS of the 
model. 

The instability of the paramagnetic solution is due to 
the presence of the quadratic term in the Hamiltonian. If 
this is missing, as for example in the spherical 3 + p spin 
glass model in which the two-body interaction is replaced 
by a three-body interaction, the paramagnetic solution is 
stable everywhere in the phase space, similarly to what 
happens for the spherical p-spin model without a field? 2 " 



B. The One Step Replica Symmetry Breaking 
Solution (1RSB) 

The stability of the RS solution q — does not depend 
on /ip. However, from the analogies with the spherical p- 
spin spin glass model we expect that for fi p large enough 
a solution with a non vanishing order parameter of the 
1RSB type might lead to a thermodynamically more fa- 
vorable phase. 

The 1RSB solution corresponds to group the n repli- 
cas into n/m clusters of m replicas. Any two replicas 
a f3 within the same cluster have overlap q%, whereas 
replicas in different clusters have overlap go < Qi ■ As 
a consequence the n x n q matrix breaks down into 
(n/m) x (n/m) blocks of dimension m x m. If the el- 
ement q a j3 with a ^ (3 belongs to one of the diagonal 
block then q a p — qi, otherwise q a p — go. The overlap 



matrix for the 1RSB ansatz can be conveniently written 
as: 

<?a/3 = (1 - qi) S a p + (qi - q ) e a/ 3 + q (23) 
where the matrix e is defined as 

1 if a and (3 are in a diagonal block 



otherwise 



(24) 



By plugging this form of q a /3 into the Eq. (jl 1|> one ob- 
tains: 



2 lim -G[q) 

n— >0 n 



5(1) ~ 9{li) + m [g(qi) - g(q )] 



x(Qo) 

TO — 1 



— lnx(<7o) 

TO 



Inx(gi). 



where, for later convenience, we have defined^ 



X(«i) 
x(Qo) 



(25) 



(26) 
(27) 



m(qi - go)- 

The saddle point equations for q and qi in the limit 
ri — * 0, obtained either from (|16f) or directly from sta- 
tionarity of (|25|l with respect to variations of go and gi, 
read 



A(g ) - 
A(gi) - 



go 



x(go) 2 
A(go) - 



= 



gi - go 
x(gi)x(go) 



o 



(28) 
(29) 



The solution of these equations depends on the value of 
to that, in the limit n — > 0, is restricted to the interval 
< m < 1. In principle any value of to which leads to a 
stable 1RSB solution can be chosen. However in the spirit 
of the saddle point calculation performed to evaluate the 
free energy we choose for any value of /i p 's the value of m 
which minimize the functional G[q] & This leads to the 
additional equation 



g(gi) -g(go) + 



i 



go 



™x(go) x(go) 2 



(gi - go) 



+ 



In 



x(gi) 



x(go) 



o. 



(30) 



The stability analysis of the 1RSB saddle shows that 
in the limit n — > the 1RSB solution is stable as long as 
the 1RSB eigenvalues^ 



A 



(i) 



A 



(3) 



-S(go) + 



1 



x(gi) 2 
l 

x(go) 2 



(31) 
(32) 



are both positive. 

The saddle point equation l|28|) admits always the solu- 
tion go = 0. It may also have solutions with < go < 1, 
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however by using arguments similar to those that lead 
to the inequality l|22[l for the RS solution, one can show 
that in absence of external field any 1RSB solution with 
<7o > is unstable since it has a negative A 3 ^ 

The 1RSB saddle point equations for q\ and m can 
be solved for any p using the same procedure used for 
the spherical p-spin spin glass model. The first step is 
to obtain g(q\) from equation (|30|) [with q a = 0] and di- 
vide it by <7i A(qi). Then using the saddle point equation 
(|29|1 [with qo = 0] to express A(qi) one ends up with the 
equation 



9iA(gi) 



= z{y) 



where 



z(y) = -2y 



l-y + lny 

(i-y) 2 



(33) 



(34) 



is the auxiliary z-function introduced by Crisanti and 
Sommers^ (CS) for the solution of the spherical p-spin 
spin glass model, and 



X(gi) 
x(Qo) 



1 - qi + m qi 



< y < 1. 



(35) 



By using y and m as free parameters eqs. I|29[) [with qg = 
0] and l|33|l can be solved for (fip,^)- A straightforward 
algebra leads to: 



/'2 



[l-z(y)] 



(p-2) qP-^l-q^l-q. + mq,) 

p (l-y + my) p 

rol 1 -^)]' (36) 



(p — 2) m 2 y{\ — y)' p 2 
1 \pz(y)-2] 



(p-2) n+mqi) 
(1 - y + my) 2 [pz(y) - 2] 



m 2 y 



(P-2) 



(37) 



By fixing the value of m in the interval [0, 1] and varying 
y these equations represent the parametric equations of 
the so called m- lines in the (/j p , \x 2 ) plane. By definition 
y can take any value between and 1 included, however, 
from Eq. I|37|) we see that since z(0) = 0, {12 becomes 
negative for y sufficiently close to 0. Setting ^2 = from 
Eq. one gets 



pz{y)-2 = {) 



(38) 



which gives the minimum value y m i n of y. The CS z- 
function (|34J) is a monotonous increasing function of y 
varying in the range < z < 1, as a consequence, \x p is 
always non-negative. 

A second condition on y comes form the stability anal- 
ysis of the 1RSB solution. A simple inspection shows 
that A^ > Aq 3 '' so that the condition which marks the 



(3) 

limit of the stability of the 1RSB solutions is Aq = 0, 
i.e., 



1 



M2 = 



(1 - qi +mqi) 2 
Using now Eq. (|3*7|) one gets the equation 
pz{y)-2-{p-2)y = 



(39) 



(40) 



whose solution gives y max , the maximum value of y for 
the 1RSB solution. 

Both boundary values y mul and y max are functions of 
p only. For example, for p = 3, we have 



while for p = 4 



0.354993.. 



fju 



0.195478. 



Vn 



1 



0.389571. 



(41) 



(42) 



The fact that for p = 3 the maximum is y max = 1 makes 
the 2 + 3 model different from any other 2+p model with 
p > 3, as we shall see in a while. 

From the stability condition l|39(l it follows that if the 
quadratic term in the Hamiltonian were missing, as for 
the already mentioned 3 + p model, then the 1RSB so- 
lution would be stable everywhere. In the Appendices ICl 
and IdI wc shall show that indeed in this case the 1RSB 
solution is the only possible non-trivial solution, besides 
the RS solution. 



1. The Transition Lines between the Paramagnet and the 
IRSB-Glass phase 

To find the transition lines which bound the 1RSB 
phase we start by noting that the m-lines do not cross and 
that the value of n 2 for which the 1RSB becomes unsta- 
ble increases as m decreases, see Eq. Q39I1 . Moreover, all 
m-lines start from \x 2 = 0. As a consequence the first m- 
line one encounters in moving from the RS phase at fixed 
fi2 < 1 and increasing fi p is the m-line with m = 1. This 
line, which marks the transition between the RS (para- 
magnetic) phase and the 1RSB (glass) phase, starts on 
the /j,p axis at the point /i p = H P (y m m, m — 1) and goes 
up to the point ji p — fi p (y max , m — 1) and = 1, as can 
be easily seen from Eq. (|39|) evaluated for m = 1. 

The transition between the RS (q — 0) and the 1RSB 
(q% 7^ 0, qo = 0) phases is not due to an instability but 
occurs because the 1RSB solution leads to a thermody- 
namically more favorable state. Since we are dealing with 
the replica trick, this means that the 1RSB solution yields 
a value of the free energy functional (|15f) larger than the 
RS solution^! 

This mechanism resembles that of ordinary first order 
transitions, and indeed the order parameter q\ jumps dis- 
continuously from zero to a finite value, and vice- versa, at 
the transition. However, the free energy remains contin- 
uous across the transition -at m = 1 the free energies of 
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FIG. 1: The RS and 1RSB phases for the 2 + 4 model in the 
(ut, H2) plane. The thick lines are the transition lines between 
different phases. 



the two solutions are equal- 
in its first derivatives. 



and no discontinuity occurs 



The 1RSB solution becomes unstable when A 



(3) 




0. 



This leads to a second transition line whose parametric 
equation in the (/Zp,/^) plane is obtained by setting y = 
2/max into eqs. lj56l - (|3*7ll and varying m from 1, to 0. For 
?7i — ► the values of both /_t p and /J2 diverge but 

M2 (1 ~ 2/max) p ~ 2 [2 - (p - 2)y max ] 



lim 



- 2) y max - p 



(43) 



and hence the 1RSB phase does not cover the full "low 
temperature" phase of the model. 

In Figure n the transition lines found so far are shown 
together with the m-lines with m — 0.7, 0.5. In the figure 
p = 4, but any p > 3 leads to a qualitatively similar 
scenario. 

The case p — 3 is special because inserting y — y m ax — 
1 into eqs. i|36|) - (|37|) one ends up with 



A*3 = m 
A*2 = 1 



< m < 1 



(44) 



Along this line q\ = 0, see e.g. Eq. and the 1RSB 

solution reduces to the RS solution. We have seen that 
the RS solution becomes unstable for /xg = 1 thus the 
critical line l|44|) marks the transition between the RS 
and the 1RSB phases. The transition is continuous in 
both the free energy and the order parameter q±. The 
transition lines for the 2 + 3 model are shown in Figure 
13 In conclusion the 2 + 3 model presents only one "low 
temperature" phase of 1RSB type and, in this respect, is 
equivalent to the spherical p-spin spin glass model in a 
field* 



C. The One-Full Replica Symmetry Broken 
Solutions (1-FRSB) 

From Figure ^ one clearly sees that for p > 3 the RS 
and 1RSB solutions do not cover the whole phase space 



FIG. 2: The phase diagram of the 2 + 3 model in the (1^3,^2) 
plane. The thick lines are the transition lines between the 
RS and 1RSB phases. In particular the thick dashed line, the 
m-line with m = 1, is the discontinuous transition while the 
horizontal thick full line is the continuous transition. 



of the 2 + p model. In the region where both the RS 
and 1RSB solutions are unstable the organization of pure 
states has a more complex structure which cannot be 
described by a simple 1RSB ansatz which groups them 
into n/m equivalent clusters. Therefore, to describe this 
region one must allow for clusters of different type. To 
this end each one of the n/m clusters is divided it into 
m/mi sub-cluster of size m\. If the procedure is repeated 
R times, dividing at each step the smallest clusters into 
yet smaller clusters, one has the i?-RSB ansatz in which 
the replica symmetry is broken i?-times. The overlap q a p 
between two replicas depends on the number of divisions 
separating the offspring clusters to which the replicas a 
and belong from the common ancestor cluster. 

A simple way proposed by Parish to parametrize the 
overlap matrix q a p for R steps in the replica symmetry 
breaking consists in dividing q into successive boxes of 
decreasing size p r , with po — n and pr+i — 1, and as- 
signing the elements q a p of the matrix q so that 

q a p = q a ni3=r = q r , r = 0,'--,R + l (45) 

with 1 = qR+i > qu . > ■ ■ ■ qi > Qo- The notation aC\(3 = r 
means that a and belong to the same box of size p r 
but to two distinct boxes of size p r +i < Pr- 

Inserting this form of q a p into Eq. one gets with 
standard manipulations 

2 R 

-G[q] = g(i) + y2(pr-Pr+i)g(q r ) 



■ln(l-q R )+J2 — In 



<lr 



t=0 



Pr q r +i 



(46) 



39,40 



where q r is the Replica Fourier Transform of q a /3t 

R+l 

Qr = Ps ^ ~ ^-i) ( 4? ) 
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The number R is arbitrary. Setting R — or R = 1 
one recovers respectively the RS and the 1RSB expres- 
sions, the latter with m = p%, while for R — > oo one gets 
the oo-RSB solution or Full Replica Symmetry Broken 
(FRSB) solution. In this limit the differences p r +i — Pr 
become infinitesimal and the set of overlaps {go, ...,qn^\ 
is replaced in the limit n — > by a non-decreasing con- 
tinuous function q{x) defined on the interval x £ [0, 1]. 

The free energy functional i|46|) for the Parisi i?-RSB 
ansatz can be conveniently expressed by using the func- 
tion 

R 

x(q)=Po + ^2(Pr+i-Pr)0(q-qr) (48) 

which equals the fraction of pair of replicas with overlap 
q a p less or equal to q. With this definition, and replacing 
the sums by integrals, one obtains, after a little of algebra 

-G[q] = f dqx(q)A(q) 

n Jq 

r qn da 

+ F77T^ +ln(1 - gfi) (49) 

Jo } q dq'x{q') 

This expression is valid for any R, and hence also for the 
FRSB solution. In the limit R — *■ oo, q r becomes contin- 
uous and we can define q(x) as the inverse of x{q). It can 
be shown that dx{q)/dq gives the probability density of 
overlaps ^Sai 

It is easy to verify that taking for q(x) = or q(x) — 
qi 9(x — to) the above functional reduces to those found 
with the RS and 1RSB ansatzs, respectively^ 

The FRSB solution with a continuous q(x) was intro- 
duced to describe the spin-glass phase of the SK model^ 
and since then it has been found in many other related 
models. A continuous order parameter function q(x) is, 
however, not general enough to describe the state of the 
2 + p model with p > 3 in the whole parameter space. 
From the stability analysis of the 1RSB solution we see 
indeed that the instability occurs because the eigenvalue 

(3) 

Ag vanishes. This eigenvalue is associated with fluctu- 
ations that involve the overlaps of one cluster as a whole 
with the other clusters as a whole^ Roughly speaking 
these fluctuations are similar to fluctuations in the RS 
phase with single replicas replaced by the clusters of to 
replicas considered as single entities. As a consequence 
we expect that, as it happens for the fluctuations in the 
RS phase, a non zero overlap go between clusters would 
stabilize the fluctuations. The solution however cannot 
be of 1RSB type since we have seen that any 1RSB solu- 
tion with go 7^ is unstable. 

Based on the analogy with the instability of RS solu- 
tion with clusters playing the role of single replicas, it 
turns out that the correct ansatz for the 2 + p model 
with p > 3 is a mixture of 1RSB and FRSB, which we 
have called 1-FRSB solution^ described by a discontinu- 
ous order parameter function in the interval [0, 1] of the 



] r- 

q(x) - 




x 



FIG. 3: Schematic form of the order parameter function q(x) 
in the 1-FRSB phase. 

form 

q (x) ^ ( 91 ^ x > m 
I q(x) tor x < m 

where q(x) is a non-decreasing continuous function in the 
semi-open interval x £ [0, to), with lim a ._ >m - q(x) = qo < 
qi and q(0) = 0^ see Fig. |31 For q = q\ one recovers 
the FRSB solution^ 3 - In appendix [D] we show that the 
1-FRSB solution is the only other possible non-trivial 
solution, beside the 1RSB (modeling a mean-field glass) 
and the FRSB (modeling a spin-glass) ones, for the 2+p 
model with p > 3. It is interesting to note that this so- 
lution also follows by solving numerically in the whole 
interval [0, 1] the Parisi equations derived from the sta- 
tionarity of the functional l|49l) with respect to order pa- 
rameter function q{x), see Appendices [XI iBl The partial 
differential equation is solved numerically by means of a 
pseudo spectral technique, see e.g. Ref. [44|, without 
fixing a priori any special ansatz for q(x). 

The free energy functional for the 1-FRSB ansatz can 
be obtained either by inserting the explicit form (15011 of 
q(x) into Eq. Q49[). or by taking in Eq. I|46|l pr = to, 
q R = qi, q R -i = q and p R - p R -\ finite as R — > oo. In 
both cases one ends up with 

2 

-G[q] = g(l)-g(qi)+m[g(qi)-g(q )] 
n 

/•So rqo f ]„ 

+ / dqx(q)A(q)+ / 
Jo Jo XKQ) 

Tfi — 1 I 

+ ln(l-gi) + -mx(<Zo) (51) 

TO TO 

where 

rqo 

X (q) = 1 - «i + m (ft - g ) + / dq' x(q') (52) 

Jq 

Stationarity of the free energy functional $ with re- 
spect to q(x) and gi leads to the 1-FRSB saddle point 
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equations: 



A(«) 



dq 

x 2 (q)' 



and 



A(gi) -A(go) 



< q < qo 



Qi - 9o 



x{Qx)x{Qo)' 



(53) 



(54) 



Finally maximization with respect to m leads to the ad- 
ditional equation 



g(qi) - s(<?o) 



1 



mx(qo) 



A(g ) 



(91 - go) 



m z 



*(?i) 



x(<7o) 



(55) 



The 1-FRSB saddle point equations fity and l|53|l are 
formally equal to the 1RSB saddle point equations 1 (29(1 
and (|30() and hence can be solved for any p with the help 
of the CS z-function. Indeed by using eqs. (|5*4*jl and ffify 

it is easy to verify that 



2 g(gi) ~ fl(go) " (gi - go) A(g ) = , . 
(gi - g ) [A(gi) - A(g )] ZU/J 

where z(y) is given by Eq. ((34(1 and y is defined as 
X(?i) _ 1 - 9i 



(56) 



V = 



x(qo) l-qi+m (qi - q ) 



< y < 1. (57) 



The saddle point equation ((53(1 is not easy to use as 
it stands. Differentiating both sides with respect to q 
to eliminate the integral one gets the more manageable 
form: 



£(<?) 



l 



x{qf 



(58) 



Equations (J^U, fB^|l and fBl^l evaluated for q = q can 
be solved for (fi p ,fj,2) as function of m and t = qo/qi- 
After a straightforward algebra one ends up with 



= 



[l-y + my(l~t)]P 
m 2 y(l - y)P~ 3 {l - t) 
1 



x 



[1- (p- l)tP- 2 + [p-2)tP- 1 ] 
[1 - y + my(l - t)} 2 
m 2 y(l — t) 2 
[y{\-tP^)(p-\)(\-t)tP- 2 ] 



(59) 



(60) 



[1 - (p - l)tP- 2 + (p - 2)V>- 1 ] 
where for any < t < 1, y is solution of the equation 

p(l - t) [1 - (p - l)tP- 2 + (p - 2)tP- x ] z(y) 

- [p - 2 - pi + p< p_1 - (p - 2)i p ] y-2+p{p- l)t p - 2 

- 2p(p-2)i p - 1 + (p-l)(p-2)< p = (61) 



Equations H59|H61|) are the parametric equations of 
the 1-FRSB m-lines which are drawn in the (fi p , pi) plane 
by fixing the value of m in the interval [0, 1] and varying 
t from to 1. The 1-FRSB to- line begins for t = at the 
boundary with the 1RSB phase and ends for t — 1 where 
the 1-FRSB solution goes over to the FRSB solution. 
The values of y for this limiting cases are: 



i = => pz(y)-2-(p-2)y = (62) 
t = 1 => y=l (63) 

By comparing Eq. Ij62(l with Eq. (|40[) one recognize that 
the value of y for t = is equal to the maximum allowable 
value of y for the 1RSB m— lines. As a consequence 1- 
FRSB TO-lines and the 1RSB TO-lines with the same m 
match continuously at the transition point between the 
two solutions. 

By evaluating Eq. H58|l for q — qo it is easy to see that 

the eigenvalue A 3 ^ [Eq. (|32|l ] is identically zero in the 
whole 1-FRSB (and FRSB) phase, in agreement with the 
marginal stability of FRSB solutions^ The eigenvalue 
A^ 1 -* [Eq. (|31|l ] remains positive in the whole 1-FRSB 
phase and vanishes for t = 1 where the 1-FRSB solutions 
disappears in favor of the FRSB solution. 

The continuous part q(x) of the order parameter func- 
tion can be obtained from Eq. (J5SJ. Indeed from this 
equation it follows that 



1-qi + to (gi - g ) + 



dq'x(q') 



VM2 + m p (p - i)q p -' 2 



(64) 



which differentiated with respect to g leads to the sought 
solution 



x(q) 



(p-l)(p-2)qP- 3 



[^ ^2 + ^, p {p-l)qp- 2 ] i,2 ' 



0<g<g O - (65) 



We note that as g — ► the probability density of the 
overlaps dx(q)/dq goes as q p ~ 4 so that it diverges for 
3 < p < 4, is finite for p = 4 and vanishes for p > 4 (see 
Fig. EI). 

Unlike the SK case^i the function q(x) for the 2 + p 
model with p > 3 is not a linear function of x for x <C 1 . 
From the solution 165(1 it is easy to see that 

q(x) ~ x 1/(p - 3 \ x->0 (66) 

so that only for p = 4 one recovers a linear behavior^ As 
a consequence dq(x)/dx vanishes for x — > for 3 < p < 4 
and diverges for p > 4. 

The function q(x) for a generic p can be obtained by 
expanding the r.h.s of Eq. ((65(1 in powers of q p ~ 2 and 
then inverting the series. As an example we give the first 
few terms for the case p = 44£ 



3/2 

q(x) = ^x 
3^4 



7/2 
^2 ^3 

6/4 



10 11/2 
13^2 5 

72/4 X 



323^2 7 
" 1296/4 ^ 



4025^ 9/2 
10368/t; 



^x 9 + O(x 10 )(67) 
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FIG. 4: The probability distributions of the overlap are plot- 
ted for different values of p. The values for the other pa- 
rameters are fi2 = 2, fi p — 3 and go = 0.2. The qualitative 
picture for small g's only depends, anyhow, on the value of 
p. Referring to the vertical ax, from top to bottom the P P (q) 
for p = 3.5, 4, 4.5, 5 and 5.5 are reproduced. At p = 3.5 (top 
curve) the P(q) diverges at q = 0, for p — 4 it goes to a finite 
value. When p > 4 it tends to zero as q — » 0. 




0.4 0.6 

FIG. 5: x c — x(qo) versus the qi — qo in the 1-FRSB phase 
for p = 4 and m = 0.5, 0.7, 1. 



and p = 5: 



£?(x) = 




(68) 



The continuous part of the order parameter functions 
ends for q — qo at the point x c = x(qo). In the 1-FRSB 
phase x c is always smaller than m and becomes equal to 
it at the boundary line with the FRSB phase. In FigureEl 
we show the value of x c as function of the difference q\— qo 
for a fixed value of m. For values of x between x c and m 
the order parameter function q(x) remains constant and 
equal to qo, and then jumps to qi as x goes through m, 
see Figure El 



q(x) 0.3 




FIG. 6: The order parameter function q(x) in the 1-FRSB 
phase of the 2 + 4 model. In the figure m = 0.8, [14 = 3 and 
= 1.4. 



1. The Transition Lines among the Amorphous Phases 
(1RSB, 1-FRSB and FRSB) 

We have seen that the 1-FRSB m-lines are the contin- 
uations into the 1-FRSB phase of the 1RSB m-lines. As 
a consequence, as the 1RSB m-line with m = 1 marks 
the transition between the 1RSB phase and RS phase, 
so the 1-FRSB m-line with m = 1 marks the transition 
between the 1-FRSB phase and the FRSB phase. The 
transition is discontinuous in the order parameter since 
<7i — qo does not vanish at the transition, but the discon- 
tinuity appears for m = 1 and the free energy and its 
derivatives remain continuous across the transition. 

The 1-FRSB m-line with m = 1 ends at the critical 
point (t = 1) 



P 



27 (p-l)(p-2)(p-3)P- 3 
.3> 



(P-2) 



where 



Qo = Qi = 



P -3 
P 



(69) 
(70) 

(71) 



For /i 2 > n\ the transition between the 1-FRSB phase 
and the FRSB takes place continuously in the order pa- 
rameter function with q\ — qo — * and x c — > m at 
the transition. The continuous transition between the 
1-FRSB and the FRSB phases occurs on the line of end 
points of the 1-FRSB m-lines. Inserting t = 1 into eqs. 
(|59|l - H6Uf) one easily gets the parametric equations of the 
critical line: 



[1-2 



2{p-3 + 3m) p 



27 m 2 [p - 
P (P - 3 H 



-l)(p- 
3 m) 2 



2)(p-3)f- 



27m 2 (p- 2) 



(72) 
(73) 
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FIG. 7: The static phase diagram of the 2 + 4 model in the 
(fi4,fJ-2) plane. 
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FIG. 8: Borders lines of the 1-FRSB and FRSB phases for 
p = 3.01 and p = 3.1. 



where < m < 1. Along this line x c = to, = and 



<7o = Qi 



P 



p — 3 + 3 m 



< to < 1 



(74) 



Finally the 1-FRSB phase is bounded by the transition 
line with the 1RSB phase. Indeed by setting t = into 
eqs. <(59l - l(6Tj)) one recovers the parametric equations of 

(3) 

the 1RSB instability line: A = and qo = 0. The tran- 
sition is continuous in both free energy and order param- 
eter function since qo — > continuously as the transition 
line is approached from the 1-FRSB side. 

All the transition lines, together with the m-lines with 
m = 0.7 and m = 0.5, and the phases of the 2+p model 
with p > 3 are shown in FigurcQ In the figure p = 4, but 
the phase diagram does not change qualitatively with the 
value of p, provided that it remains larger than 3. 

In the limit p -> 3 the 1-FRSB and FRSB phases 
shrink to zero while the transition lines separating the 
two phases collapse smoothly onto the vertical line (0, 1x2) 
with fi2 > 1 and the horizontal line (1, ji p ) with < fi p < 
1 where go = <Zi = 0, see Figure [S] One then smoothly 
recovers the phase diagram of the 2 + 3 model, Figure |21 

From Figure |S] we see that the continuous transition 
line between that 1-FRSB and the FRSB phases displays 
a point of vertical slope in the (fXp,^) plane. Along 
the continuous transition line between the 1-FRSB and 
FRSB phases the point of vertical slope is attained for 



m 



(00) 



2 p- 3 

3 p- 2 



where 



(00) 



(p-2) (Ap-3)P 



6{p 
p(p- 



1) ( P -3)p- 1 

2) (4p - 3) 2 



12 ( P - 3) 2 



(75) 

(76) 
(77) 



and qo = q± = (p — 2) /p. This point exists for any p > 3. 



Similarly the point of infinite slope along the transition 
line between the 1-FRSB and the 1RSB is attained for 



(00) _ ^ (1 — Umax) 

{p - 2) y max 



(78) 



where j/ max is given by the solution of Eq. I|40|) . For this 
value one has 



(oo) = (p-2) (4p-3)P 

Mp 6(p- vip-sy- 1 ' 

(00) ^max(l 1/max) 



l<2 



4(p- 2)f" 2 



(79) 
(80) 



and qo = 0, 91 = (p — 2)/p. This point exists only for 
3 < p < 3.5197.... 



D. The Full Replica Symmetry Broken Solution 
(FRSB) 

For the FRSB solution the order parameter function 
q(x) is continuous. The equations for the FRSB phase 
are easily obtained from those of the 1-FRSB by setting 
go = Qi and m = 1 so that only the continuous part 
of the order parameter function survives. In the FRSB 
phase the function x(q) is still given by Eq. H65|l but with 
qo = Qi solution of [see Eq. 1 (55)1 ] 



S(gi) = 



1 



(l-qi) 2 



(81) 



The order parameter function q(x) in the FRSB is 
shown in Figure |5J 

By defining r = — 1 from eqs. I|81[) and 1 (65(1 it follows 
that when the RS instability line \ii = 1 is approached 
from the FRSB side then 



'7i 



T 

2' 



0^ 



(82) 
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FIG. 9: Order parameter function q(x) in the FRSB phase. 
In the figure for p = 4 /14 = 2 and [i% — 1.5. 



and 



(p- l)(p- 2) /TNP-3 



H 



(83) 



It is easy to see that for a generic p the expansion of qo 
in powers of r coincide with that of 1 — 1/yT+T up to 
order (3(t p ~ 2 ) not included. For example for p = 4 one 
has 



?i~^ + ^4-1)t 2 +0(t 3 ) 



(84) 



while for p = 5 



fli ~ r 

yl 2 8 



-(5 + 4 M5 )r 3 + 0(r 4 ) (85) 



and so on. 

The transition between the FRSB and RS phases oc- 
curs for \i% = 1 where both q\ and x c = ie(<Zi) vanish and 
the FRSB solution goes over the RS solution q = 0. The 
transition line (ji p , 1) ends at the crossing point with the 
m-line with m = 1. The transition between the FRSB 
and the RS phases is continuous in the order parameter 
function and hence in the free energy. 



III. COMPLEXITY 

The 1RSB and 1-FRSB ansatz both contain the pa- 
rameter m which gives the location of the discontinuity 
in the order parameter function. Strictly speaking the 
replica calculation does not give a rule to fix it. Going 
back to the expression of the moments of the replica par- 
tition function, Eq. I|10|> , one indeed sees that the replica 
calculation requires that for any overlap matrix q a p the 
free energy functional must be extremized for N — > 00 
with respect to the elements of the matrix. However, it 
does not say anything about the structure of the matrix 
q a p. This means that in the 1RSB [and 1-FRSB] ansatz 
the free energy functional $ must be extremized with re- 
spect to qi [and q(x)] but not necessarily with respect to 



to, since it is related to the matrix structure. This rises 
the question of which value of m has to be taken when 
there exist different values of to, all of which leading to 
a stable solution. In the solution discussed so far the 
value of to yielding the maximum of the free energy 3 ^ 
was chosen. 

The free energy functional, Eq. (|15f) . evaluated on the 
stable saddle point solution gives the free energy of a 
single pure state. As a consequence, choosing for to the 
value which maximizes the free energy functional is ther- 
modynamically correct, provided that the logarithm of 
the number of different pure states with the same free 
energy, called complexity or configurational entropy, is 
not extensive. If the configurational entropy is exten- 
sive, it gives a contribution to the thermodynamic free 
energy which must be considered when computing the 
extrema. In other words if the number of states is exten- 
sive the extrema of the thermodynamic free energy fol- 
low from a balance between the single state free energy 
and the configurational entropy. This is what happens in 
systems with a 1RSB phase, as first noted in the p-spin 
modelfSiiS and changes the condition for fixing the value 

of TO. 

We shall not give here the details of the direct calcu- 
lation of the complexity for the 2 + p model, but rather 
we shall use the shortcut of deriving it from a Legen- 
dre transform of the replica free energy functional with 
respect to to. 

To be more specific the complexity Sit in the 1RSB and 
1-FRSB phases is obtained as the Legendre transform of 
/3m$(m) where $(to) is the replica free energy functional 
(|15fl evaluated with the 1RSB or 1-FRSB ansatz keeping 
to as a free parameter: 



Eit(/) = max[/?m/ - (3m$(m)] 



(86) 



We shall use for the complexity the notation En to stress 
that it is obtained from the Legendre transform, and to 
distinguish it from the "self-energy" function E(x) used 
in section [H] Strictly speaking this is the complexity 
density, even if it is customary to call it just complexity. 
In the Legendre transform, Eq. (|86[) . / is the variable 
conjugated to to 



/ 



<9to<I>(to) 
dm 



(87) 



and its value equals the value of the free energy inside a 
single pure state for the given value of to. Introducing 
this expression into the Legendre transform one gets the 
following relation 



Eit(/)=/?m 



2 9$ (to) 



where to(/) is the value of to found by solving Eq. 1)87(1 . 

By using the expression l|51|) for the functional G[q] the 
complexity of the 1-FRSB solutions of the 2 + p model 
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reads: 



2E lt (m) 



-to 2 [g(q t ) - g(q )] 
Qi - go 



In 



x(qo 



Kilo) 



m (gi - go) A(g ) 



where go and Qi must be evaluated as function of /Lt p , 
and to using the saddle point equations (|54|l and (|58|1 . 
Alternatively we can use Eq. I|54ll to eliminate to in favor 
of gi so that the expression of the complexity for the 1- 
FRSB solutions becomes: 



2Eh(?i) = 1-ln 



[A(gi) - A(go)](l - gi) 2 
(<?i - go) 
(gi - go) 

1 A(gx)-A(g ) 
- U - Qi) 



l - gi gi - go 

.9(91) - 5(90) - (gi - go)A(gi) 



[A( gi )-A(g )] 2 
where go is given by the solution of 

(1-gx) 2 [A(g 1 )-A(g )] : 



M2 + M P (p- I)?? 



(gi - go) 5 



(89) 



(90) 



and gi is such that ^(gi) = /, i.e., m(qi) = m(f). 

The complexity for the 1RSB solution is obtained just 
setting go — into the 1-FRSB complexity [and neglect- 
ing Eq. H9Ufl ]. A simple check of the 1RSB complexity 
consists in verifying that for [12 = one recovers the 
complexity of the spherical p-spin models 

By varying qi one selects 1RSB or 1-FRSB solutions 
with different to. As a consequence not all values of gi 
between and 1 are allowed but only those which lead 
to stable solutions must be considered. This means non- 



negative eigenvalues A^ and Aq 01 for 1RSB solutions and 

non-negative eigenvalue A^ for 1-FRSB solutions. The 

(3) 

eigenvalue A is identically zero for 1-FRSB solutions. 
The requirement that only solutions with non-negative 
A^ are physically acceptable is also know as the Ple- 
fka's criterion42iSi Here it comes out naturally from the 
stability analysis of the replica saddle point, however it 
can be shown to have a more general validity. 

The complexity Sit is the logarithm of the number of 
states of given free energy, divided by the system size N. 
It is, therefore, clear that in the thermodynamic limit 
only solutions with a non-negative complexity must be 
considered. All others will be exponentially depressed 
and hence are irrelevant. 

The static solution discussed in previous Sections was 
obtained by imposing d<&(m)/dm = 0. The complexity 
is consequently zero for the static solution, and the 
number of ground states is not extensive. 

The solution with the largest complexity, of both 1RSB 
or 1-FRSB type, is the one for which A^ vanishes, 
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FIG. 10: Eit(gi) as function of q\ in the region where only 
1RSB solutions are stable and have non-negative complexity. 
Only the part of the curve Ei t (qi) > is shown. The thick line 
shows the relevant part of the curve. The dashed line corre- 
spond to unstable solutions. The point marked "Dynamic" is 
where A^ vanishes and corresponds to the solution of largest 
complexity. The figure is for the 2 + 4 model with /12 = 3 and 

fl 4 = 10. 
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In Figures HU1 1111 and 1 1 21 we show the behavior of Si t (gi) 
in the three relevant regions where (i) only 1RSB solu- 
tions have non-negative complexity and are stable, (ii) 
both 1RSB and 1-FRSB solutions have non-negative 
complexity and are stable and (iii) both 1RSB and 1- 
FRSB solutions have non-negative complexity but only 
the latter are stable. 

The condition of maximal complexity (|91|) is known as 

the "marginal condition" since for A^ 1} = the saddle 
point is marginally stable. 

In the relaxation dynamics the eigenvalue A^ is re- 
lated to the decay of the two-times correlation function 
to the "intermediate" value gi, and hence the marginal 
condition comes naturally in as the condition for critical 
decayi 2 ^ For this reason the solution of maximal com- 
plexity is also called the "dynamic solution" as opposed 
to the "static solution" discussed so far which, on the 
contrary, has vanishing complexity. 

In the FRSB phase A^ is identically zero and the two 
solutions, static and dynamic, coincide. 



IV. THE DYNAMIC SOLUTION 

In this paper we shall not give here the full derivation 
of the dynamic solution, and of the marginal condition 
l|9ip. starting from the relaxation dynamic equations but 
rather we shall rely on the fact that the dynamic solu- 
tion can be obtained from the replica calculation just 



13 



E,,^) 0.004 



V 3) >o 



Dynamic 




(1.65 



0.75 



FIG. 11: En(gi) as function of q\ in the region where both 
1RSB and 1-FRSB solutions are stable and have non-negative 
complexity. We have indicated explicitly which part of the 
curve corresponds to each solutions. Only the part of the 
curve Ei t (gi) > is shown. The thick line shows the relevant 
part of the curve. The dashed line correspond to unstable so- 
lutions. The point marked "Dynamic" is where vanishes 



and corresponds to the solution of largest complexity, 
figure is for the 2 + 4 model with fi2 — 3 and [14 = 7. 



The 



dynamic solutions. 



A. The Dynamic 1RSB Solution 

The equations of the dynamic 1RSB solution are given 
by Eq. (|29[1 with qo = and by the marginal condition 
(|91|l . Solving these equations for {p p ,H2) and using qi 
as a free parameter one gets the parametric equations of 
the dynamical 1RSB m-lines 



M2 



(p-2)qr 3 (l-qi) 2 (l-qi+mqi) 

(p- 2)(1 - gi) - mgi 
(p-2)(l-q 1 )2(l-q 1 + mq 1 ) 



(92) 
(93) 



For any < m < 1 the maximum allowable value of q\ is 
fixed by the requirement that [ii > 0, while the minimum 

by the requirement that the eigenvalue A3 , Eq. P2*|l 
with q = 0, be non-negative. In the dynamic approach 
this eigenvalue controls the long time relaxation of the 
two-times correlation function 2 ^ and hence must be non- 
negative. A straightforward calculation shows that 
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FIG. 12: Eit(gi) as function of qi in the region where 1RSB 
and 1-FRSB solutions have non-negative complexity but only 
the latter are stable. We have indicated explicitly which part 
of the curve corresponds to each solutions. Only the part of 
the curve Ei t (gi) > is shown. The thick line shows the 
relevant part of the curve. The dashed line correspond to 
unstable solutions. The point marked "Dynamic" is where 
Aj 1 ' vanishes and corresponds to the solution of largest com- 



plexity. 
fl 4 = 4. 



The figure is for the 2 + 4 model with /i2 = 3 and 



using the marginal condition instead of stationarity of 
the replica free energy functional with respect to mw^^ 
It can be shown that this shortcut applies also to the 
1-FRSB solution^ 

The static and dynamic solutions differ only for what 
concern the 1RSB and 1-FRSB phases, therefore here 
we shall only discuss shortly the main differences in the 
phase diagram which follows from the 1RSB and 1-FRSB 



p-3 

with < m < 1. 
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P 



p-2 



(94) 



1. The Dynamic Transition Line Between the Paramagnet 
and the IRSB-Glass Phase 

The dynamic transition line between the RS and the 
1RSB phases is given by the dynamic 1RSB m-line with 
m = 1. Inserting m = 1 into eqs. I|92|l - (|93|l one obtains 
the parametric equations of the transition line 



1 



M2 = 



(p-2)<$- a (l- qi ) 
(p-2)-(p-l)gi 
(p-2) (l-gija 



(95) 
(96) 



with (p_3)/(p-2)<gi < (p-2)/(p-l). In the (fi p ,» 2 ) 
plane the line begins on the fx p axis at the point 



(P-1) 



p-i 



(p - 2)p- 2 ' 
and goes up till the point 

(p-2r 2 







Mp 



(p - 3)P" 3 : 



^2 = !• 



(97) 



(98) 



where Ag '' vanishes. The transition between the RS and 
the dynamic 1RSB phase is discontinuous in the order 
parameter q\ since it jumps from zero, on the RS side, 
to a finite value on the m-line with m = 1. 
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The dynamic 1RSB phase is bounded by the critical 



line of equation A 



(0) 



which marks the transition 



between the dynamic 1RSB and the dynamic 1-FRSB 
phases. The explicit form of the equation of this transi- 
tion line is obtained by inserting qi = (p— 3)/(p— 3 + to), 
see Eq. lj93)l. into the equations of the dynamic 1RSB m- 
line and reads 



/'2 



(p - 3 + m)P 



(p - 2) 2 (p - 3)p- 3 m 2 
(p - 3 + m) 2 



(p — 2) 2 TO 2 



(99) 
(100) 



with < to < 1. 

As expected, the dynamic transition lines do not coin- 
cide with the static ones but, in the (/ip,^) plane, they 
are displaced toward lower values of fi p with respect to 
the corresponding static transition lines, see Figure IT51 



B. The Dynamic 1-FRSB Solution 

The equations of the dynamic 1-FRSB solution are 
given by the saddle point equations (|5*5|) and ijoTI) and 
by the marginal condition 191|) . As a consequence, the 
parametric equations of the dynamic 1-FRSB TO-lines are 
still (|59fl - l|tj0|) but with the value of y given by 



l-jp-^tP- 1 + {p-2)tP~ 
p-2-(p- Vjt + tP- 1 



(101) 



where < t = qi/qo < 1. The continuous part of the 
order parameter function is given by Eq. (|65() with [see 
Eq. 



qa=tqi 



t(i-y) 



1 — y + my (1 — f) 



(102) 



The dynamic 1-FRSB m-line are drawn in the (/j, p , ^2) 
plane by fixing the value of to and varying t from to 1 . 



1. The Dynamic Transition Line between the 1RSB and the 
1-FRSB Amorphous Phases 

By setting t — into the equations of the dynamic 
1-FRSB TO-lines and varying to from 1 to one recovers 
the critical line qo = and A3 = which marks the 
boundary with the dynamic 1RSB phase. Indeed for t = 
Eq. ((TUTjl yields y = l/(p - 2) so that eqs. 
(|60|l reduce to the parametric equations l|99|) - (|100|l of the 
critical line. Moreover on this line qi — (p— 3)/(p— 3+to), 
the same value found from the dynamic 1RSB solution, 
therefore as it happens for the static solution the dynamic 
1-FRSB TO-lines and the dynamic 1RSB m-lines with the 
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same to match continuously on the critical line A 
[and go = 0], This transition is continuous since qo goes 
to zero as the transition line is approached from the 1- 
FRSB side while q\ is continuous through the line. 
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FIG. 13: Phase diagram /42-P4. Counterclockwise the RS, 
1RSB, 1-FRSB and FRSB phases are plotted, separated by 
the static phase transition lines (full curves). The PM/1RSB 
and the 1RSB/1-FRSB transitions also occur in the dynamics 
and the relative lines are drawn as dashed curves. Their con- 
tinuation as FRSB/1-FRSB transition lines are the dynamic 
(dashed) and static (full) m = 1-lines, computed in the 1- 
FRSB ansatz. They merge at the "End Point" (see inset). 
For a comparison, we also plot the dynamic and static m- 
lines with m = 0.5. They merge on the FRSB/1-FRSB phase 
transition line above the end point. As m decreases from one, 
the whole continuous FRSB/1-FRSB line is covered. Inset: 
the discontinuous transitions FRSB/1-FRSB (^12 > 1) and 
PM/1RSB (jja < 1). 



The dynamic 1-FRSB m-line with m = 1, continua- 
tion of the 1RSB m-line with m = 1 into the 1-FRSB 
phase, marks the boundary between the 1-FRSB and 
FRSB phases. Along this line the order parameter is dis- 
continuous since q\ jumps from zero in the FRSB phase 
to a non-zero value on the line. The discontinuity oc- 
curs at m — 1 so that the free energy remains continuous 
despite the jump in the order parameter. The dynamic 
1-FRSB m-line with m = 1 starts from the end point 
(|98ll of the dynamic 1RSB m-line with m = 1 and stops 
at the same end point (|69[) - H70I) of the static 1-FRSB m- 
line with m = 1. From this point the transition between 
the 1-FRSB phase and the FRSB phase occurs continu- 
ously in the order parameter, i.e., with q\ — qo - * as 
the transition line is approached form the 1-FRSB side 
(see Fig. |T3|). 

The continuous transition between the 1-FRSB and 
FRSB phases occurs along the critical line obtained by 
setting t = 1 into eqs. (|59|I - H60I) and varying m from 1 to 
0. From Eq. ((TUT)) it follows that 



y = l- P —^{l-t) + 0{{l-tf), 



t 



r 



(103) 



so that the end point of the dynamic 1-FRSB m-line co- 
incides with the end points of the static 1-FRSB m-line 
for any to, and not only for m = 1. Therefore, the dy- 
namic and the static continuous critical lines between the 
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FIG. 14: Phase diagram T-J4 in J2 units. Clockwise a param- 
agnetic phase (PM) , a glassy phase (G-I) described by means 
of a 1RSB ansatz, another glassy phase represented by a 1- 
FRSB ansatz and a purely spin-glass phase (SG) computed in 
the FRSB ansatz occur. The transitions to and from the two 
glass-like phases take place both as dynamic (dotted lines) 
and as static (full lines) . In temperature the dynamic transi- 
tion lines are always above the relative static lines. 



1-FRSB and the FRSB phases coincide. Indeed, when 
studying the static solution we have seen that along the 
continuous transition line between the 1-FRSB and the 
FRSB solutions the eigenvalue A^ vanishes so that the 
difference between the two solutions disappears. On this 
line both solutions have zero complexity and it remains 
equal to zero in the whole FRSB phase. 

In Figure ED we show the full phase diagram of the 
spherical 2 + p spin glass model in the (/z p ,/i2) space 
with both the static and the dynamic critical lines. 

By noticing that both fi p and /12 are proportional to 1 
we see that the discontinuous dynamic transition occurs 
at a temperature higher than that of the equivalent static 
transition, as can be clearly seen from Figure 1141 where 
the phase diagram in the T/J2 = and J p /J2 = 

•\/2/ip/ ' {pp2) plane is shown. 



V. CONCLUSIONS 

In this paper we have provided a detailed study of the 
phase space of the spherical 2 + p spin glass model us- 
ing the static approach of Ref. |2] which employs the 
replica method to evaluate the disorder-averaged loga- 
rithm of the partition function. By performing the Leg- 
endre transform of the replica free energy functional we 
have defined the complexity function that, whenever it 
is extensive, counts the number of equivalent different 
metastable states. This allowed us to discuss the dynamic 
solution as the solution which maximizes the complexity. 
In both solutions, i.e., static and dynamic, the model 
displays four different phases, characterized by different 



replica symmetry breaking schemes, in which the system 
can find itself as the thermodynamic parameters and the 
intensity of the interaction is changed. One of the nice 
feature of this model is that it can be completely solved 
even in the phase described by a full replica symmetry 
breaking. To our knowledge this is the first example of 
an analytical FRSB solution. 

The main result of the present study is the existence 
of two phases exhibiting the qualitative features of glassy 
materials. One is the 1RSB solution, displayed by those 
systems that are considered as valid mean-field models 
for the glass state. Even though no connection with the 
microscopic constituents of a real glass former can be 
set, the collection of spins interacting through a p-body 
quenched disordered bound behaves very much like, e.g., 
the set of Si02 molecules in a window glass. The single 
breaking occurring in the 1RSB solution corresponds, in a 
dynamic interpretation£2iSi£4 to a time-scale bifurcation 
between the fast processes in a real amorphous material 
(the /3-processes) and the slow processes (a) responsible 
for the structural relaxation. 

Another phase emerges in the study of the spherical 
2 + p spin glass model. Something not occurring in any 
Ising spin glass models In a whole region of the phase 
space the stable phase is, indeed, described by means 
of an overlap function q(x) that is continuous up to a 
certain value qo and then displays a step, as in aforemen- 
tioned 1RSB solution. We call it the one step-full RSB 
solution (1FRSB). Exploiting the static-dynamic analogy 
once again, in this phase, in the relaxation towards equi- 
librium, a first time-scales bifurcation takes place ( "a — j3 
bifurcation") just as above, but it is no more unique. As 
the time goes by a continuous set of further bifurcations 
starts to occur between slow and even slower processes, 
as in the case of a proper spin-glass^ In the continu- 
ous part, any kind of similarity between (ultrametrically 
organized) states is allowed but above go the hierarchy 
ends in only one extra possible value: the self-overlap, or 
Edwards- Anderson order parameter q\. 

The stability of the 1-FRSB solution in the replica 
space is not limited to a single, self-consistent, choice of 
the order parameter (in particular the point discontinuity 
can change in a certain interval) and this implies that in 
each point of the phase diagram belonging to this phase 
there will be an extensive number of metastable states, 
having free energies higher than the equilibrium free en- 
ergy. In order for this phase to appear a strong enough 
couple interaction must be present (that is the source 
of the continuous, spin-glass like contribution) but the 
p-body interaction must have a broader distribution of 
intensities than the 2-body (J p > J2). 

If this new phase can be considered as a glassy phase 
different from the 1RSB one, the phase diagram that we 
have computed describes an amorphous-amorphous tran- 
sition, with the second glass having a much more compli- 
cated structure outside the single valleys (for very long 
time-scales in the dynamic language). Whether there is 
a correspondence with the amorphous-amorphous tran- 
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sition between hard-core repulsive and attractive glassy 
colloidsi^i^iSiiSiiS and at which level the analogy can 
be set is yet to be clarified and further investigation of 
the dynamic properties is needed in order to make a link 
between the interplay of 2-body and p-body interaction 
in spherical spin glasses and the role of repulsive and at- 
tractive potentials in colloids undergoing kinetic arrest. 



APPENDIX A: THE PARISI EQUATION AND 
THE SOMMERS-DUPONT FORMALISM 



The Parisi equation for the 2 + p spin glass model is 
more easily obtained starting from the free energy func- 
tional in the replica space written as 

Pf[q, A] = -s(oo) - lim — 9{<l a f}) 

a/3 

+ ! irn o2^E A ^^ (Al) 

a/3 

- lim - lnTrv cxp [ 1 V \ a /3 o a u^ 
V «V J 
where the matrix X a p is the Lagrange multiplier associ- 
ated with the replica overlap matrix q a p, see Eq. <|14[) . In 
particular the diagonal element X aa — A is the Lagrange 
multiplier that enforces the spherical constraint q aa = 1 . 
Stationary of f[q,X] with respect to variations of X a p 
and q a/ 3 leads to the self-consistent equations (a ^ (3) 

X a p = A(q a p) (A2) 

Tr CT o a o fi exp (j2 a p cr a a 



Tr a exp \ J2 a p X a /3 cr a a>- 



(A3) 



By applying the Parisi's replica symmetry breaking 
scheme an infinite number of times and introducing the 
functions X(x) and q(x), < x < 1, one for each matrix, 
the free energy functional (|A1|) for the spherical model 
can be written as 

Pf = - s (oo)-^[g(l)-J dxg(q(x))] 

+ \\X-f Q dxX^x)]-^-^-. 



+ OC 



dy 



exp 



2A(0) 



</»(0,y) (A4) 



where cf>(0, y) is the solution evaluated for x = of the 
Parisi equation 



<j>{%, y) = K x ) W'{ x , y) + x <f>'( x , y) 2 ] 

with the boundary condition 

1 v 2 



2 A(l)-A 



(A5) 



(A6) 



In writing the Parisi equation we have used the standard 
notation in which a dot "0" denotes the derivative with 
respect to x while the prime "</>"' the derivative with 
respect to y. 

The advantage of this equation is that it can be solved 
numerically without specifying a- priori the form of q(x) 
and X(x). The first problem one is facing when solving 
the Parisi equation is that the functional (|A1|I must be 
extremized over all possible solutions of the Parisi equa- 
tions, which can be numerically uncomfortable. This, 
however, can be overcome using the Sommers-Dupont 
formalism. 56 The idea is to introduce a different func- 
tional whose value at the stationary point coincides with 
the extrema of the free energy functional (|A1|) over all 
possible solutions of the Parisi equations. 

This is easily achieved by introducing the Parisi equa- 
tion into the functional via the Lagrange multiplier 
P(x,y). The new functional is hence 



Pfv = Pf + 



+ 00 



dx 



dyP(l,y) 

+ 00 



0(1, y) 



l 



y 



2A(1)-A 



dyP(x,y) (j>{x,y) 



X(x) </>"(x,y)+x<f>'{x,yf 



(A7) 



where (3f is the free energy functional (|A1|> . 

The functional f3f v is stationary with respect to vari- 
ations of P(x,y), P(l,y), <j)(x,y), <fi(0,y), the order pa- 
rameter functions q(x) and X(x) and A. Stationarity with 
respect to P{x,y) and P(l,y) just gives back equations 
(|A5I) and (| A6|> . while stationarity with respect to 4>(x,y) 
and 0(0, y) leads to the differential equation for P(x, y): 

P(x, y) = ^X(x) [P"(x, y) ~ 2x [P(x, y) m(x, y)] ' 

(A8) 

where m(x,y) = <fi'(x,y), with the boundary condition 



P(o,y) 



exp 



y2^A(0) V 2A(0) 
S(y), as A(0) 



(A9) 



It can be shown that P(x, y) is the probability distri- 
bution of the local field y at the scale q(x). 

Finally stationarity with respect to q(x), X(x) and A 
gives 



X(x) = A(q(x)) 



(A10) 



q(x) 



+ 00 



dyP{x,y)m{x lV ) 2 (All) 



and 



A(l)-A 



1 



\-q{l) 



(A12) 
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FIG. 15: The numerical order parameter function q(x) for the 
2 + 4 model with /12 = 2 and /14 = 3. The 1-FRSB solution is 
clearly seen. 



P(x,y) 0.4 




FIG. 16: Field probability distribution p(x, y) for the 2 + 4 
model with y>2 = 2 and /i4 = 3 evaluated from the numerical 
resolution of the Parisi-Sommers-Dupont equations at x = 0.4 
and x = 0.5. The underlying full curve (barely visible) is the 
analytic solution given by Eqs. <B9L IIB141 . 



From Eq. (|A1 l|l and the identification of P(x,y) with 
the local field distribution it follows that m(x, y) repre- 
sents the local magnetization at scale q(x) in presence of 
the local field y. It obeys the differential equation 



m(x,y) = -^( x ) i m "( x >y) + 2xm(x,y)m'(x,y)} 



with initial condition 



1^(1, y) = [1 - q(l)]y 



(A13) 



(A14) 



The partial differential equations l|A8(l and (|A13I) can 
be solved numerica lly usi ng the pseudo-spectral method 
developed in Refs. |44l46| . In Figure IT51 we show the or- 
der parameter function q(x) found solving the equations 
for the 2 + 4 model with ^ 2 = 2 and fi 4 = 3. The 1-FRSB 
structure is clearly seen. In Fig. 1161 both the numerical 
and the analytical solutions for P(x, y) are displayed. 



APPENDIX B: SOLUTION OF THE 
PARISI-SOMMERS-DUPONT EQUATIONS 

It is not difficult to write the Parisi-Sommers-Dupont 
partial differential equations for the 1-FRSB case, cfr. 
Ref. With obvious notation, the Parisi-Sommers- 

Dupont functional with the 1-FRSB ansatz reads 



Pfv = Pf + 



+00 



dyP(m,y) (f>(m,y) 
,,2 



2 Ai — A — m (Ai — Ao) ■ 

/+00 _ 
dyP(x,y) 4>(x,y) 
-OO 

h(x) U"(x,y)+xcf>'(x,y) 2 ]] (Bl) 



where [3f is the free energy functional 



Pf = s{oo) 

1 



[g(l)-(l-m)g( qi ) 



dxg(q(x))] 



1 - f m 

^ (! ^ m ) - / dx\(x)q(x) 
z Jo 



2 



■J-m 

2m 



A x - A 



Ai — A — m (Ai 
- dy 



exp 



Ao) 



2A(0) J 



0(0, y) (B2) 



Stationarity of (|B lfl with respect to variations of 
P(x, y) and P(m, x) yields the Parisi equation (|A5|) with 
x restricted to the interval [0, m] and initial condition 



<p(m,y) 



1 



2 Ai — A — m (Ai — Ao) 



(B3) 



Similarly stationarity with respect to variations of 
4>(x, y), 0(0, y), q(x), X(x) for < x < m gives again eqs. 
!|A8|). (|A^|) . I|A10() and (|A11|I . while stationarity with re- 
spect to A leads back to eq (| AI 2|> with A(l) replaced by 
Ai. 

For the 1-FRSB ansatz there are three more equations. 
The first two follow from stationarity with respect to q\ , 
Ai and read 



Ai = Afa) 



Ai — Ao 



<?i - <7o 



(1 - qi)[l - qi - m{qi - q )} 



(B4) 



(B5) 



Finally stationarity with respect to m, the discontinu- 
ity point in the order parameter function, leads again to 
equation (|55|) . It is not difficult to recognize in these 
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equations the saddle point equations for q$, q\ and m 
derived for the 1-FRSB phase. 

For the spherical model the Parisi-Sommers-Dupont 
equations can be solved analytically. Indeed defining 



F(x) =Ai- A-m(Ai- A ) - / dx' x 1 \{x') (B6) 

J X 

it is easy to verify that the solution reads 
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F(x) 



i(x,y) 



dx 



F{x') 



F(x) 



and 



P(x,y) 



y / 2ira(x) 2 



exp 



where 



a(xf = F{xf 



dx 



2a(xY 



, A(x') 



(B7) 



(B8) 



(B9) 



F(x') 2 

Inserting Eqs. (|B8|1 - (|B9|) into Eq. (jAllf) one gets 



(BIO) 



q(x) 



dx' 



F(x'f 



(Bll) 



Taking the derivative of eqs. I|B6|) and IjBllf) with respect 
to x it is easy to show that 



(B12) 



which integrated yields the equality x( x ) = ^/F{ x )- 
Thus Eq. IjBllfl becomes 



q{x) 



dx' X(x') x{x') 



'\2 



(B13) 



Inverting this relation one gets back Eq. (J53J. We note 
that changing the integration variable from x to q in Eq. 
(|B11|) one can show that F 2 (q) = £(q) so that the vari- 
ance a(x) 2 of the local field distribution 1B9I) can be writ- 
ten as 



a(xf = T,(q(x))q{x) 



(B14) 



This completes the solution of the Parisi-Sommers- 
Dupont equations for the spherical 2+p spin glass model 
in the 1-FRSB phase. The FRSB solution can be ob- 
tained taking go — ► Qi and m — > 1. 



APPENDIX C: THE SPHERICAL s+p SPIN 
GLASS MODEL 



The spherical s+p spin glass model is the generaliza- 
tion of 2 + p models in which the two-spin interaction is 
replaced by a s-body spin interaction: 



UN 



UN 



ii<...<i s 



n<...<ij, 



.( C1 ) 

In the following we shall assume that s < p, even if the 
Hamiltonian is trivially invariant under the exchange of 
s and p. 

The study of the s+p model follows closely that of the 
2+p with the replacement of g(x) by [see Eq. (|12fl ] 



/ \ s ftp p 

g(x) = — x H — -x y . 
s p 



(C2) 



As it happens for the 2+p models the RS solution with 
q =/= is unstable, however, at difference with these, 
the RS solution with q = is stable everywhere in the 
(fip, ft s ) plane since for q = the relevant eigenvalue Ai, 
ea. (|21|l . is identically equal to one for any s > 2. 

For fi p and/or fi2 large enough a thermodynamically 
more favorable 1RSB solution appears. This solution has 
qo = and q\ given by the saddle point equation 129|) 
[with #0 = 0]. For the static solution the value of m is 
fixed by Eq. l(3Tj|) . or equivalently by Eq. (j23J), which fol- 
lows from stationarity of the free energy functional with 
respect to variations of m. These equations can be solved 
for any s and p with the help of the CS z-function and one 
obtains the parametric equations of the static m-lines: 



ftp 



fts 



p (1 - y + my) p [2 - sz(y)] 



(p — s) m 2 y(\ — y)P 2 2 

s (1 - y + my) s [pz(y) - 2] 
(p — s) m 2 y(\ — y) s ~ 2 2 



(C3) 
(C4) 



where y = (1 — <Zi)/(l — q\ + mqi). Notice that, as 
expected, the expressions are symmetric under the ex- 
change of s and p. 

The 1RSB solution never becomes unstable since for 

("3") ____ 

s > 2 the eigenvalue A , Eq. (|3T)l , evaluated at go = 

is always positive, and the eigenvalue A^ , Eq. J22Jl, 
remains positive for any ft p > 0. 

As a consequence the limit on y are given by the con- 
ditions 



ft P = o z (y» P ) = - 

2 

ft s =Q => z(y^ s ) = - 



(C5) 
(C6) 



The CS ^-function is an increasing function of y therefore 
if s < p then y^ p < y fJ-3 . For s = 3 and p = 4 we have 

= 0.195478..., ^ 3 = 0.354993... (C7) 
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FIG. 17: The phase diagram of the 3 + 4 model in the (/j.3, (14) 
plane. The thick dashed lines are the m-line with m = 1, i.e. 
represent the discontinuous transition between the RS and 
the 1RSB phases. 



To find the dynamic transition line between the 1RSB 
and the RS phase equation (|33[) must be replaced by the 
marginal condition Jjnj. Using qi as the independent 
variable it is easy to derive the parametric equations for 
the dynamic m-lines: 



1 (a - 2 + m) gi - (s - 2) 

1 (p - 2) - (p - 2 + m) gi 

(P-*) gr 2 (l-9i) 2 (l-9i + ^9i) 



(C8) 
(C9) 



The range of gi is fixed by the requirement that both \x v 
and fi s be non-negative. This yields the boundary values 



Hp = =>■ qi = 



s-2 



Hs = =>■ gi = 



-2 + 771 

p-2 



P- 



(C10) 
(Cll) 



In Figure 1171 we show the phase diagram for the s + p 
model with s > 2 and p > s. In the figure s = 3 and 
p = 4, however any choice of s > 2 and p > s leads to a 
qualitatively similar phase diagram. 



APPENDIX D: THE RSB SOLUTIONS 

In this Appendix we show that the spherical 2 + p spin 
glass model admits only solutions of 1RSB, FRSB or 1- 
FRSB type. Theprocedure is the same as that of Ap- 
pendix 2 of Ref. 0. 

We start from the free energy functional for the -R-RSB 
ansatz [see Eq. (|46|) ] 



-G[q] 

n 



dq x(q)A(q) 
dq 



J^dq'x(q') 



+ ln(l-« R ) (Dl) 



where 



x(q) =Po +22(p r +i -p r )0(q- <h)- (D2) 



The saddle point equations are obtained by varying the 
above functional with respect to x{q): 



SG[q] = f dqF{q)5x{q) (D3) 
Jo 



where 



F(q) = A(q) 



dq' 



tdq"x{q") 



(D4) 



and 



Sx i(l) = ^2 ~ S Pr) d (l " Qr) 

R 

-^2{Pr+i-Pr)6(q-q r )6q r (D5) 



r=0 



By requiring stationarity of G[q] with respect to the q r 
and the p r one gets, respectively 



F(q r ) = 0, 



dqF(q) = 0, 



= i2 



R 



(D6) 



(D7) 



The function F(q) is continuous, thus Eq. (|D7|I implies 
that between any two successive q r there must be at least 
two extrema of F(q). If we denote these by q* , then the 
extremal condition F'(q*) = implies that 



dqx(q) = — , 



(D8) 



The left hand side of this equation is a concave function, 
i.e. with a negative second derivative, since x(q) is not 
decreasing with q. For the 2+p model E(g) = /j 2 +/J p (p— 
l)q p ~ 2 so that the right hand side is concave for small q, 
provided p > 3, and convex for large q, see Figure 1181 
For p = 3 the right hand side is convex. 

As a consequence Eq. (|D8|) admits at most two solu- 
tions and hence only one step of replica symmetry break- 
ing is possible. Moreover it must be go = 0. Indeed, from 
Eq. (|D6|> and from the fact that in absence of external 
field F(0) = it follows that 



F(0) = F(q ) = F( qi ) = 



(D9) 



which would imply the presence of at least three extrema, 
which is not possible. 

Up to this point the conclusions do not differ much 
from those found for the p-spin model. Here however the 
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FIG. 18: Typical behavior of l/-^/E(g) as function of q for 
p > 3. 

presence of a concave part in the right hand side of Eq. 
(|D8(1 for p > 3 makes possible different solutions. 

Equations (|D6(I and (|D7|I can be solved by a continuous 
replica symmetry breaking solution with F(q) — in a 
given range of q since in this case both equations would 
be identically valid. For this solution Eq. (|D8(I must 
also be identically valid and this can only be true for 



< q < qo , where go is the value of q for which 1 / y/Y.(q) 
changes concavity, where both sides of the equation are 
concave function of q. Indeed if Eq. 1|D8|) were valid also 
for values of q where the right hand side is convex this 
would imply that x(q) is a decreasing function of q, which 
is not possible since dx{q) / dq is the probability density of 
the overlaps. For the same reason for p = 3 a continuous 
replica symmetry breaking solution is not allowed since 
in this case the right hand side of Eq. (|D8|I is purely 
convex. The same applies to the s+p models with s > 2, 
so that also these models admits only a 1RSB phase. 

The possibility of a continuous replica symmetry 
breaking solution for q < qo does not rule out the pres- 
ence of discrete replica breakings with q r > qo. These 
additional discrete breakings must satisfy eqs. I|D6(I and 
(|D7(1 . Therefore since the right hand side of Eq. ljD8|) 
is convex for q > qo arguments similar to those which 
leaded to the conclusion that for the 2 + 3 model only 
1RSB solutions are possible show that at most only one 
more discrete break with qi > qo is possible. 

No other possible non trivial solutions exist for the 
spherical 2 + p spin glass model so we conclude that the 
model admits only solutions of 1RSB, 1-FRSB or FRSB 
type. 
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